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We consider charge relaxation in the mesoscopic equivalent of an RC circuit. For a single-channel, 
spin-polarized contact, self-consistent scattering theory predicts a universal charge relaxation resis- 
tance equal to half a resistance quantum independent of the transmission properties of the contact. 
This prediction is in good agreement with recent experimental results. We use a tunneling Hamil- 
tonian formalism and show in Hartree-Fock approximation, that at zero temperature the charge 
relaxation resistance is universal even in the presence of Coulomb blockade effects. We explore 
departures from universality as a function of temperature and magnetic field. 
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There is increasing interest in the dynamics of meso- 
scopic structures motivated by the desire to manipulate 
and measure quantum phenomena as rapidly as possi- 
ble. It is thus of great importance to characterize the 
time scales governing the electron dynamics in simple 
mesoscopic structures. An elementary but fundamental 
building block is the quantum coherent capacitor Q . As 
in the classical case, the low frequency dynamics of a 
mesoscopic capacitor is determined by a charge relax- 
ation time Tfic- For a quantum coherent capacitor the 
RC time can still be written as the product of a resis- 
tance and a capacitance, i.e. trc = RqCf^. However, 
due to the coherent nature of electron transport through 
mesoscopic structures, both the electrochemical capaci- 
tance C^, which determines the charge on the capacitor 
and the charge relaxation resistance Rq, which governs 
the charge fluctuations, now crucially depend on coher- 
ence properties of the system. The capacitance is 
related to the imaginary part of the AC conductance but 
it can also be obtained by the differentiation of a thermo- 
dynamic (grand-canonical) potential [1, H, Q IE Hi - The 
charge relaxation resistance Rq is related to the real part 
of the AC conductance and therefore requires a dynamic 
theory. In analogy to the classical RC circuit depicted in 
the upper left corner of Fig. ^ one has for the mesoscopic 
system 



(1) 



This equation will be taken as a definition of and 
Rq . If the cavity-reservoir connection permits the trans- 
mission of only a single spin polarized channel, a self- 
consistent scattering matrix approach, gives at zero tem- 
perature a resistance equal to half a resistance quan- 
tum Q 



R - ^ 



(2) 



We emphasize that the factor of 2 is not connected to spin 
but is rather due to the fact that the cavity connects 
to only one electron reservoir. More astonishing, even 



counter-intuitive, is the fact that Eq. (2) is independent 
of the transmission properties of the channel. 

In a seminal experiment, J. Gabelli et al. have re- 
cently measured both the in and out of phase parts of 
the AC conductance of a mesoscopic RC circuit. In their 
experiment one "plate" of the capacitor consists of a sub- 
micromcter quantum dot (QD) and the other is formed 
by a macroscopic top gate. The role of the resistor is 
played by a tunable quantum point contact (QPC) con- 
necting the QD to an electron reservoir. The results of 
this experiment are in good agreement with the theo- 
retical predictions of 0, in particular they confirm the 
universality of the single channel charge relaxation resis- 
tance by using a strong magnetic field to spin polarize 
the electrons. However, it is a priori unclear whether the 
results derived in Ref. yl still hold in the presence of sin- 
gle charge effects 0, 0,0 which must become important 
if the transmission through the QPC becomes small. In- 
deed the experiment observes Coulomb blockade oscilla- 
tions of the capacitance as a function of the gate voltage. 
It is the aim of the present work to present a theoreti- 
cal description for the charge relaxation resistance in the 
presence of Coulomb blockade effects. 

The mesoscopic RC circuit along with the principal 




FIG. 1: Schematics of the mesoscopic capacitor. A cavity is 
connected via one lead to an electron reservoir at voltage V{t) 
and capacitively coupled to a backgate with voltage Vg. The 
coupling matrix elements Tmn are defined in the text. The 
inset shows the corresponding classical RC circuit. 
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model parameters is shown in Fig. V{t) = V^c cos((jjt) 
is the time dependent voltage applied to the electron 
reservoir while Vg is the voltage applied to the gate and 
C is the geometrical capacitance between the QD and 
the gate. The matrix elements Tmn, to be defined below, 
describe the indirect coupling between the dot states m 
and n via the reservoir. 

Our first goal is to show that in the single channel 
case Rq is universal also in the Coulomb blockade regime. 
For that purpose we treat the cavity at the Hartree-Fock 
(HF) level [10. The starting point of our calculation is 
the relation 



/(0 = -e|3{tr[G<(t,t)]}, 



(3) 



which expresses the tunneling current through the QPC 
as a function of the lesser (matrix) Green function (GF) 
G^(t,t') of the dot. We are interested in the regime 
where the modulation energy is small compared to the 
level spacing in the dot. Treating the dot as zero dimen- 
sional, the Hamiltonian of this system is 



H = Ht 



Hr 
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+ h.c.), (4) 



= (l/27r) / dEc-'^^^-*'^ f{E), is the Fourier transform of 
the Fermi function. = 2TTpLt'^*t'^ are the couphng 

matrix elements in the wide band limit |l2l |. Here and 
in the following, we use the matrix notation A^a.na' = 
^mn^<r<T' i which takcs advantage of the fact that spin is 
conserved in Q). An important property of the coupling 
matrix elements is that ^mn^ki — ^mi^km from which it 
immediately follows that for arbitrary matrices A and B 



tiiV' AT" B] = tr[r'"yl]tr[r'"B], 



(7) 



where the trace is over a basis of dot states with spin 
a. Since we are interested in the linear conductance, we 
expand © to linear order in V and U and find after 
double Fourier transformation 

G<{E, E') = G%{E)T<{E, E')Gfg{E') + 0{U^V'), 

(8) 

where rf{E,E') is the double Fourier transform of 
T<{t,t') = iT{l + icj){t,t'))f{t - t') with 0(t,t') = 
4>v{t,t') - (j)u{t,t'). With © and © the linear re- 
sponse tunneling current at frequency to becomes I{uj) — 
g{uj)[V{Lu) - U{uj)], with 



where Hl - Eka Ekait)cl^Cka, with Eka{t) - + 

eV{t), describes the (non-interacting) electrons in the 
isolated lead and fj. ^ is the tunneling matrix clement be- 
tween a reservoir momentum state k and the n-th single 
particle dot state, both with spin a. A change in the gate 
voltage is modeled as a shift of the Fermi energy Ep in 
the reservoir and we set Vg = 0. The Hamiltonian of the 
dot reads ilu = J2na '^nadi„dna+Ec{Ndot+J^{t))^ ■ Here 



Ec = e^/2G is the electrostatic charging energy, Ndot = 
J2ma '^ma'^rncr is the particlc number in the dot and 
eN{t) — GU{t) gives the polarization charges between 



dot and gate produced by the time dependent voltage at 
the reservoir This polarization charge in turn, leads 
to a time-dependent (Hartree) potential U{t) inside the 
dot and we may write Hd = Y.na ^ncr{t)d'l^dncr+EcN^^t, 
with imcrit) = tma + eU{t). In HF approximation, the 
retarded (advanced) GF of the dot takes the form pLall^ 



G«(^)(t,t') =e''^^(*'*')GfJ^)(i-i'), (5) 

where Geq^\t-t') = G^^^\t,t')\v^^^o is the equilibrium 
retarded (advanced) HF Green function and (j)ij{t,t') ~ 

dTlKj). The equal time lesser GF is obtained via the 
Keldysh equation 

G<(t,t)= \ dtx f dt2G'^{t,ti)T<{ti,t2)G^ih,t). (6) 



Here r<{t,t') = iFe^"^^ <*'*')/(*- <') with 0y(i,i') = 
/j* dTV{T) is the lesser coupling self-energy and f{t — t') 



/ dEFiE,L,MGfgiE)TGtq{E-u;)], (9) 



where i^(i?, w) = [/(_E+ijj) — /(£')]/a' and we have set ?i = 
1. To obtain the AC conductance G(w) = I{uj)/V{uj), we 
need the internal potential U{iu). For this we note that 
in the present single lead system, the displacement cur- 
rent —iujeN{w) is equal to the tunneling current so that 
g{ijj)\V{ijj) — U{uj)] = —iujCU((jj) and consequently |0| 
U{uj) = g{uj)V {lo) / [—iuC + git-o)]. Expanding the con- 
ductance to second order in frequency, we then obtain 
after restoring the units 



R„ = 



h JdEfjEMDjE)^] 
2e2 {JdEf'{E)tr[D{E)]y 



(10) 



where /' = df/dE and D{E) = Gf^{E)TGfg{E). Using 
0, we have tT[D'^{Ef] = tT[D''{E)]'^ and hence at zero 
temperature 



Ra 



(E^MEf)? 



(11) 



where Vcr{E) = tr[£''^(i?)]/27r is the density of spin a 
states in the dot. Eqs. H9I11() are central results of this 
work. In particular, Eq. demonstrates that for a 
single (spin polarized) channel Rq is still given by Eq. |5J). 

In the following, we use Eq. to investigate the 
magnetic field dependence of Rq. We consider here a 
dot with two spin degenerate levels with bare energies 
eio- and £20- = eicr + A respectively. In the numerical 
calculations, we set Ed IS. = 2.5. We are interested 
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FIG. 2: Magnetic field dependence of tlie charge relaxation 
resistance Rq. The upper panels show Rq as a function of 
the Zeeman splitting Ab and the Fermi energy Ep for weak 
■K'y/Ec — 0.13 and strong ■K'y/Ec = 1.4 coupling. All energies 
are given in units of the bare level spacing A. The lower 
panels show the corresponding total dot charge. 



in the regime of low magnetic field, where the Zeeman 
splitting As = fJ-BSB < A. For simplicity, we fur- 
ther assume that rj„„ = ^mn = 7^ f^i' m-,"- £ {li2}. 
The equilibrium HF retarded GF of the dot may be 
written as G^^iE) = [G§{E)-'^ - Sg^ + iT/2]-^, with 
the non-interacting equilibrium GF of the isolated dot 
{G^{E)U^ ^ S^„{E - e„, + iO+)-i and the HF self- 
energy 



ncr 



(12) 
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The most important feature of this self-energy is that it 
correctly excludes the unphysical self-interaction terms 
(m = n ~ I and a = a') of the Hartree approxima- 
tion and consequently leads to the appearance of the 
Coulomb gap across Ep, which is the essential spectral 
signature of the Coulomb blockade effect. The "mean 
fields" (4 

solving the set of equations 



nuydna) are determined self-consistently by 



{dlnadncr) 



1 

2^ 



dEf{E)[G^JE)rGtJE)r^„. (13) 



Because of the interaction, Vcr{E) depends on the level 
occupancies and we must distinguish two cases. Solving 
Eq. H13|) numerically in the strongly blockaded regime 
TT^/Ec ^ 1, we find that Rq is non-universal even as 
B ^ (upper left panel of Fig. 13). This is due to the fact 
that the dot charge is strongly quantized in this regime, 
as shown in the lower left panel of Fig. El which leads 
to a gap of order Ec + A^crcr' + (1 — 5cya')^B between 
the highest occupied state with spin a and the lowest 



unoccupied state with spin a' . There are thus four well 
separated (separation Ec), narrow (width ~ 7) reso- 
nances in the total density of states as a function 
of Ep. We can understand the particular behavior of Rq 
shown in the upper left panel of Fig.|21at specific values 
of the Fermi energy. If Ep is on a resonance with spin 
CT, then v„i{Ep) = i'cr{Ep)5aa' ■ Therefore on resonance, 
i.e. for Ep e {ei^., e2a, ei?? + £^c + As, e2CT + -Bc + As}, we 
expect to have Rq = /i/2e^ according to Hll() . Further- 
more, in the middle between two consecutive resonances 
Va{Ep) = i^(Ep) and so Rq takes on its minimal two 
channel value h/Ae^. In the opposite limit of strong cou- 
pling 1 < TT^/Ec, the Coulomb blockade is smeared out 
and the charge on the dot is not strongly quantized any- 
more (see lower right panel in Fig. \^ . Thus at very low 
magnetic fields, the two spin states are nearly simulta- 
neously charged and the degeneracy is not appreciably 
lifted. Therefore, Rq = h/Ae^ at low magnetic fields as 
shown in the upper right panel of Fig. [3 As A^ in- 
creases the broad resonances in the density of states cor- 
responding to the two spin degenerate levels split into 
consecutively overlapping resonances and Rq starts to 
increases except in the middle between two neighboring 
resonances. Finally, a crossing of the two innermost reso- 
nances is observed for As = 0.75. Such a crossing occurs 
when the Zeeman splitting approaches the effective level 
spacing /S. + Ec{{nia) — {n^a)), which interestingly is seen 
to be smaller than the bare level spacing here. This effect 
is peculiar to the strongly coupled regime and is due to 
enhanced exchange interactions, which favor the conse- 
quent population of states with equal spins IHIii. At 
the crossing point the densities of both spin states are 
equal and Rq takes on its minimal value. 

As a second application, we investigate the tempera- 
ture dependence of Rq in the high magnetic field limit, 
where the incoming electrons are effectively spin po- 
larized and there is only a single transmitting channel 
through the QPC. We consider here a quantum dot with 
two (bare) levels ei and £2 = ci + A and suppress the now 
supcrfiuous spin index. In the low temperature regime 
kgT ^ 7, we may expand v{E) and v{E)'^ around Ep 
assuming that these functions vary slowly in the range 
Ep ± kBT/2. This yields, to first non- vanishing order in 




ly'iEf 



v{Ef) 



(14) 



where v'{E) = Thus in HF approximation, the 

lowest order correction is proportional to the square of 
the energy derivative of the density of states at the Fermi 
energy. This explains the presence of the peaks seen in 
Fig. El to the left and right of the two resonances at 
Ep = ei = 1 and Ep ^ t\ = ei + i?c + A, for the two 
lowest temperatures (i — 100 and 13 = 12.5. The fact that 
for /? = 12.5, Rq docs not identically vanish at resonance 
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FIG. 3: Charge relaxation resistance Rq as a function of 
Fermi energy Ef for different temperatures. The lower three 
curves are for a two level dot. The upper most curve is the 
high temperature asymptote. The inset shows Rq as a funtion 
of temperature for Ef — t\. 



where v'[Ep) = 0, is due to higher order terms in the 
low temperature expansion, which involve non-vanishing 
higher order derivatives of v. The dotted horizontal line 
at Rq — h/2e'^ marks the zero temperature result. In the 
opposite limit of very high temperature /cgT 3> A + i^c, 
where ~f'(E) = /3/4 and in the weakly coupled regime 
7 <C A, where the density of states is well approximated 
by a sum of displaced Lorentzians of width 7. we can 
estimate the remaining integrals in the numerator and 
denominator of ^ a.s - J dEf'{E)v{Ef = /3Af/47r7 
and -/ dEf'{E)u{E) = /3M/4 respectively, where M is 
the number of density of states peaks under the broad 
curve f'{E). Therefore at very high temperature we find 



Rn 



e2 TTjMP 



(15) 



This high temperature limit is shown as a dashed hor- 
izontal black line in Fig. O for the two highest tem- 
peratures. In the present two level system, M = 2 
and as expected the agreement with the numerical in- 
tegration of (|10|l is good for the highest temperature 
f3~^ = 50 ^ A + Ec ~ 3.5. For the intermediate temper- 
ature A = 1 < = 2 < A + i?c = 3.5, there is already 
a significant deviation from the asymptotic result. In the 
inset of Fig. 01 we show Rq as a function of the tempera- 
ture on the first resonance at Ep = ei. The dashed line 
corresponds to the high temperature asymptote H15|l and 
the dotted line is the low temperature limit 114|) . 

In this work we have analyzed the charge relaxation 
of a mesoscopic capacitor in the linear regime of coher- 
ent dynamical transport. We have shown that the sin- 
gle channel zero temperature charge relaxation resistance 
Rq is universal even in the presence of single charge ef- 
fects, described in the Hartree-Fock approximation. This 



shows in particular that charge relaxation of a quantum 
coherent capacitor is faster than one could naively expect 
based on classical arguments. We obtain the magnetic 
field dependence of Rq in the two channel case (electrons 
with spin), where we identify two qualitatively different 
regimes of weak and strong coupling. In the former, the 
degeneracy of both spin states is lifted by the interac- 
tion at all field strengths and Rq is non-universal. In the 
latter regime, the degeneracy is lifted only at finite field 
and at zero field Rq is universal and equal to its mini- 
mum two channel value h/Ae^. The finite temperature 
behavior of Rq for a two level spin polarized system is, 
to lowest order, determined by the logarithmic deriva- 
tive of the density of states with respect to energy. In 
the multilevel case the HF approximation gives a reason- 
able qualitative picture of the underlying physics. The 
important case of B = 0, for a single strongly coupled 
level in the dot, requires a treatment of Kondo physics 
and will be discussed elsewhere . 

Our work demonstrates that mesoscopic charge relax- 
ation is a physically very interesting process and provides 
a basis for the understanding of experimental data in the 
low and high magnetic field ranges. 
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